Unstable regimes for a Bose-Einstein condensate in an optical lattice 
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We report on the experimental characterization of energetic and dynamical instability, two mech- 
anisms responsible for the breakdown of Bloch waves in a Bose-Einstein condensate interacting with 
a ID optical lattice. A clear separation of these two regimes is obtained performing measurements 
at different temperatures of the atomic sample. The timescales of the two processes have been 
determined by measuring the losses induced in the condensate. A simple phenomenological model 
is introduced for energetic instability while a full comparison is made between the experiment and 
the 3D Gross-Pitaevskii theory that accounts for dynamical instability. 
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I. INTRODUCTION 

The interest in the system made by neutral atoms in 
optical lattices has constantly been growing in the last 
decades since the development of efficient laser cooling 
techniques, which opened the possibility of observing 
quantum effects on the motion of atomic ensembles. The 
physics of quantum particles in periodic potentials can be 
described in terms of Bloch waves [H and indeed many 
effects originally predicted for electrons in a lattice of ions 
have been observed for ultracold thermal atoms moving 
in optical lattices 2] and, more recently, in quantum de- 
generate samples with the observation of long lived Bloch 
oscillations in a degenerate Fermi gas 3]. In particular, 
the achievement of Bose-Einstein condensation (BEC) in 
dilute atomic gases has allowed the possibility to repeat 
these experiments with ensembles of particles all occupy- 
ing the same Bloch state, in principle enhancing the visi- 
bility of these quantum effects |j,l3,la,l3. However, when 
the density of the sample increases, as in the case of a 
trapped condensate, interactions among the atoms form- 
ing the BEC may significantly change the simple single- 
particle picture. The interaction-induced nonlinearity is 
responsible for the observation of many other phenom- 
ena, notably the phase transition from a superfluid to a 
Mott-insulator 8] and the generation of bright gap soli- 
tons 0. Furthermore, it can be shown that there exists 
a range of parameters for which non-linearity makes the 
Bloch-like solutions of the wave equation describing the 
system unstable. 

In this paper we report on the experimental charac- 
terization of the unstable regimes for a BEC in a ID 
optical lattice, obtained through a clear measurement of 
the timescales describing the evolution of the system, in 
remarkable agreement with the theory. 

The phenomenon of instability of a BEC in a ID op- 
tical lattice has already been the subject of experimen- 
M M M El and theoretical [ll El El El El 
[T^ I20I I21L |22^ works. In particular, in the framework 
of the Gross-Pitaevskii (GP) theory in a periodic poten- 
tial, there are different mechanisms responsible for the 
breakdown of the initial superfluid state. On one hand. 



as demonstrated by Landau in the context of superfluid 
helium |2^, there is a critical velocity (related to the 
sound velocity of the system) beyond which the system 
can lower its energy by emitting phonon-like excitations 
which deplete the original state. This kind of instabil- 
ity, that is closely related to the energy spectrum of the 
system, is called energetic instability and it has been ob- 
served for a harmonically trapped condensate in [23 |. 
On the other hand, due to the interplay between non- 
linearity and periodicity in the GP equation governing 
the dynamics of the system, for certain values of lattice 
height and velocity an arbitrarily small fluctuation of the 
original state may grow exponentially in time, thus de- 
stroying the initial Bloch state. This kind of instability, 
common to many non-linear systems in a periodic poten- 
tial, is usually called dynamical or modulational instabil- 
ity since it is connected to the dynamic equation which 
describes the system. In Fig. Q] we show a schematic dia- 
gram of the stability of the first band of Bloch states for 
a condensate in a ID optical lattice as a function of the 
lattice height. The theory predicts the existence of three 
regions depending on the the condensate quasimomen- 
tum corresponding to a regime where Bloch waves are 
stable solutions, and energetically or energetically and 
dynamically unstable. 

Although the distinction of these two mechanisms is 
straightforward from the theoretical point of view, in the 
experiments it is much more difficult to separate ener- 
getic and dynamical instability since, no matter which 
is the mechanism responsible for the onset of instabil- 
ity, the original BEC superfluidity will be compromised 
and losses of atoms in the ground state are expected. 
Furthermore, as it is shown in Fig.Q the stability of the 
system depends on the height of the optical lattice: when 
the lattice height is greater than the chemical potential 
of the BEC, the so-called tight binding regime (lattice 
height > 5 Erec in Fig. energetic instability is only 
a narrow boundary between the regions of stability and 
dynamical instability. On the other hand, when the lat- 
tice is weak compared to the chemical potential, as in 
the experiments reported in this paper, there is a wide 
range of values of quasimomentum for which the system 
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FIG. 1: Schematic plot of the stabihty of Bloch matter waves 
in the first band of a ID optical lattice: white area corre- 
sponds to the stable region, in the light gray area the conden- 
sate is energetically unstable and in the dar k gray area it is 
both energetically and dynamically unstable l2l| . Vertical 
axis is the Bloch quasimomentum q and horizontal axis is the 
lattice height in recoil energy (see Sec. HTj) . 



is energetically but not dynamically unstable. For these 
reasons the experimental technique used in the v ery first 
experiments devoted to the study of this regime 
namely the excitation of dipolar oscillations in a peri- 
odic plus harmonic potential, cannot be used to separate 
the effects of energetic and dynamical instability. This 
happens because, during a single experimental run, the 
condensate velocity (and therefore its quasimomentum) 
evolves throughout the first band and thus explores both 
regimes, as indicated in Fig.^ The phenomena observed 
in were attributed to energetic instability because 
the ID Gross-Pitaevskii theory used to analyze the data 
could not reproduce some of the observed experimental 
features. This in turn suggested that finite temperature 
effects could play a role. This interpretation lead to a 
debate 10] and finally in 21], by making a comparison 
with the full 3D theory, it was shown that the density 
profiles observed in [Tfl could be attributed to dynami- 
cal instability. 

In the present paper we show that the precise control 
of the BEG quasimomentum, obtained using an optical 
lattice moving at constant velocity, is crucial to distin- 
guish the different unstable regimes. The experimental 
procedure implemented in this work allows a compari- 
son with the predictions of the theory and a study of the 
stability of the condensate also in the excited bands of 
the periodic potential 13]. The structure of the paper 
is the following: in Sec. ^ we describe the experimen- 
tal set-up, in Sec. IIIII we present novel results obtained 
for energetic instability and discuss the role of thermal 
excitations, then in Sec. |^ we present a comprehensive 
investigation of the phenomenon of dynamical instability 
and show that it is possible to clearly distinguish it from 
energetic instability. In the Appendix □ we discuss the ef- 
fect of the harmonic trapping potential on the dynamics 



of the BEG adiabatically loaded in the moving optical 
lattice. 



II. EXPERIMENTAL SET-UP AND 
PROCEDURE 

We produce a BEG of ^^Rb atoms in the hyper- 
fine state \F = l;mF = —I) using a standard double 
MOT apparatus. Our condensate typically contains 
3 X 10^ atoms and it is produced by radiofrequency- 
induced evaporation of the atomic sample confined in an 
elongated magnetic harmonic trap characterized by an 
axial and radial frequency of ujz = x (8.74 ± 0.03) Hz 
and u;± = 27r x (85 ± 1) Hz respectively. The ID optical 
lattice is formed by two counterpropagating laser beams 
obtained from a Ti:Sa source and it is aligned along the 
axis of the magnetic trap. The interference profile has 
a spatial period of A/2, where A 820 nm is the wave- 
length of the laser. Using two single mode fibers we ob- 
tain two gaussian beams with a radius of 200 /im, much 
larger than the condensate radial size. In this paper we 
report measurements carried out with a lattice height of 
s = 0.2 and s = 1.15, where s is the height of the optical 
potential in recoil energies {Erec = h? / {2\^m)^ m being 
the mass of a Rb atom). The lattice height is calibrated 
via Bragg scattering of the condensate and is monitored 
throughout the experiment 25]. The frequencies of the 
two beams are controlled independently by two acousto- 
optic modulators which use the same timebase. This al- 
lows us to precisely control the velocity of the lattice vl^ 
which is related to the frequency difference between the 
two beams, Az/, by the following relation 



Vl 



A 



(1) 



In the laboratory frame the lattice potential can be writ- 
ten as 



Vl{z) = sErecCOS^ 



2 f 27r{z — VLt) 
A 



while the confining magnetic potential is 



1 



Vho{^) = ^m{u;lr^^u;lz^). 



(2) 



(3) 



After producing the condensate in the harmonic po- 
tential, we load it adiabatically into a single Bloch state 
ramping the intensity of the lattice from zero to its final 
value in a time tu :^ 2 ms. During the ramp time the lat- 
tice velocity is kept constant and no acceleration is used. 
This process can be viewed as a deformation of the spec- 
trum of the system from a free particle one (i.e parabolic) 
to a Bloch one, in which the states are labeled by a band 
index n and a quasimomentum q (that has the periodicity 
of the reciprocal lattice, 2qB = 47r/A = 2mvB/K)' If the 
ramping time is sufficiently long so that this deformation 
can be considered adiabatic, the condensate is transferred 
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to a Bloch state with quasimomentum q = \vL\m/h^ 
belonging to the first band if \vl\ < vb^ to the sec- 
ond band if vb < \vl\ < "^vb and to the n^^ band if 
{n — 1)vb < \vl\ < nvB ^4, 5] (further details are given 
in the Appendix). To observe the effect of non-linearity, 
we maintain a high density in the sample keeping the 
harmonic trapping potential on during the whole exper- 
imental procedure. This is crucial to enter the regime 
where the effects of dynamical instability are observable. 

In order to study the two mechanisms of instability 
we perform a time resolved analysis of the losses in the 
condensate induced by the lattice for different values of 
quasimomentum, band index and lattice height. The 
general procedure is the following: after loading the con- 
densate into a single Bloch state we let it interact with 
the moving lattice for a time t, then we switch off adia- 
batically the lattice, release the atoms from the magnetic 
trap and measure the number of atoms remaining in the 
condensate by taking an absorption image along the ra- 
dial direction after 28 ms of expansion. This allows us 
to reconstruct the evolution of the number of atoms as a 
function of time and to measure the lifetime of the con- 
densate r as a function of q. Without the optical lattice 
the EEC lifetime is limited to r = (23 ± 3) s by the heat- 
ing due to current noise in the coils producing the mag- 
netic trap. The introduction of a stationary lattice (i.e. 
a standing wave) introduces two further heating mecha- 
nisms: resonant photon scattering and lattice vibrations 
due to mechanical noise on the optics setting the path of 
the two beams. Using the Ti:Sa laser at a wavelength of 
820 nm (74THz detuned from the Dl line) makes heat- 
ing from resonant photon scattering completely negligi- 
ble, while we measured that the second effect becomes 
dominant when the lattice height is higher than one re- 
coil energy (r = (17 ± 2) s with s = 1.15). In the case 
of a moving optical lattice this lifetime can be strongly 
reduced if instability mechanisms are activated. 

In particular we found that, also for a very small ve- 
locity, the lifetime measurement is critically affected by 
the presence of a residual cloud of noncondensed atoms 
surrounding the BEC. Even a barely detectable fraction 
of thermal atoms can significantly shorten the lifetime of 
the condensate and thus the temperature of the system 
must be taken into account. For this reason, in order 
to control the temperature of the sample, the radio fre- 
quency used in evaporative cooling is kept on after the 
production of the BEC (RF-shield). 



III. THE ROLE OF THE THERMAL 
FRACTION: ENERGETIC INSTABILITY 

The effect of a thermal fraction can be seen in Fig. [3 
There we show two series of pictures taken at different 
lattice quasimomenta ranging from to 0.2 qB and for 
two values of the final radiofrequency used for the evap- 
oration, corresponding to different temperature of the 
atomic cloud. When we have an almost pure condensate 



(thermal fraction < 20%, limited by our imaging sensitiv- 
ity), the number of atoms remaining after 15 s of interac- 
tion is not sensitive to the lattice velocity (bottom part of 
Fig.I2|), while in case of a mixed cloud (thermal fraction 

35%) even a small velocity leads to a strong reduction 
in the number of atoms (top part of Fig. [21). Note that 
the velocities presented in these pictures are well below 
the threshold for the onset of dynamical instability for 
the experimental parameters of Fig. [3 corresponding to 
q 0.5 qB (see Sec. II V|) : we will come back to this point 
at the end of this section. This strong reduction in the 
number of atoms, triggered by the dissipation provided 
by the thermal cloud, can be qualitatively interpreted as 
the effect of energetic instability. 

In order to test this hypothesis quantitatively we have 
measured the number of atoms in the condensed fraction 
as a function of quasimomentum for three different times 
of BEC-lattice interaction, using a 35% thermal fraction. 
The results are shown in Fig. [31 As one can see, the 
number of atoms N slowly reduces with increasing q up 
to a critical quasimomentum qo^ after which it remains 
constant. The smooth behavior of A/" as a function of q 
below this critical value seems not to be compatible with 
the threshold process expected in |1^]. However this be- 
havior can be understood with an argument similar to 

IT" 

the one of reference 10]. As we already pointed out, 
the onset of energetic instability is related to the sound 
velocity within the condensate: once the center of mass 
velocity exceeds this velocity the system can lower its en- 
ergy by emitting phonon-like excitations. In the presence 
of axial confinement, as in the experiments reported here, 
the system is inhomogeneous and thus the sound veloc- 
ity, which depends on the density, is not constant along 
the direction of the lattice, vanishing at the edge of the 
condensate. One can therefore expect that, once the cen- 
ter of mass velocity v is non-zero, there exists a fraction 
of the condensate in which the local sound velocity Vs 
is lower than v and therefore that fraction becomes en- 
ergetically unstable. The sound velocity for an infinite 
cylindrical condensate in a ID optical lattice is given by 

m 

Vs = {^^^o, (4) 

where g is an effective interaction constant which takes 
into account the presence of the lattice, m* is the effec- 
tive mass for the condensate Bloch state and no is the 
peak density of the sample. Assuming that a similar re- 
lation holds for the local density in an inhomogeneous 
cylindrical condensate, we can argue that the condensate 
is locally energetically stable if Vs{z) > Vn,q where Vn,q 
is the Bloch velocity of the condensate loaded in band n 
with quasimomentum q. As one can verify, for the values 
of quasimomentum and lattice height involved in these 
measurements, the Bloch dispersion is only slightly differ- 
ent from the free particle one for which Vn^q = hq/m and 
it is correct to assume g = g = 47rfi^ a^/m and m* = m 
where m is the mass of a Rubidium atom and is the 
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FIG. 2: Absorption images of the condensate interacting for t = 15s with a lattice with s = 0.2 for different values of 
quasimomentum ranging from to 0.20 and for respectively a condensed fraction of about 65% (top) and no detectable 
thermal component (bottom). 



s-wave scattering length. The above condition for ener- 
getic stability can be recasted as 



ni{z) > 



g m 



(5) 



where we indicate with ni{z) the density of the conden- 
sate along its axis. One can therefore obtain the fraction 
of the condensate which is energetically stable by inte- 
grating the density over the region satisfying Eq. © 



fqo{Q) = ^ I I I dxdyn 

^^0 Jni>Cq^J J 



(r) 



(6) 



where A^o is the number of atoms in the condensate and 
n(r) its density distribution. 

In our experiments the lattice produces a weak mod- 
ification of the density profile. It is therefore correct to 
assume as local density (i.e. averaged over the lattice 
spacing) the Thomas-Fermi profile of the condensate in 
the magnetic potential 



n{x,y,z) = no 1 



z 



(7) 



which gives the fraction of the condensate which is ener- 
getically stable 



fqoiQ) 



where 6{x) is the Heaviside function and qo 




{gmno)/ (2Ti^) is the threshold value for energetic in- 
stability in a homogeneous cylindrical condensate with 
peak density no, namely the quasimomentum for which 
Eq. (0) cannot be satisfied for any z. The value of qo can 
be theoretically calculated also beyond the approxima- 
tion of free particle dispersion and is independent from 
the BEC-lattice interaction time provided that the den- 
sity is not significantly reduced by the losses induced by 
the instability. 

In order to derive a simple expression for A/" as a func- 
tion of q and t one has to make some assumptions on 



the decay induced by energetic instability. We assume 
that for a given q the number of atoms in the stable 
fraction (i.e. Nq fq^^q)) is constant in time, while the 
number of atoms initially in the unstable fraction (i.e. 
A'o (1 — fqoiq))) decays with a time behavior b{t). We 
thus obtain for A^ the following expression 

N{q, t) = No UM) + m ^0 (1 - Uo (9))- (9) 

Note that Nq b{t) can also be viewed as the number of 
atoms remaining in the condensate after a time once it 
is entirely unstable {fq^ = 0) and A" no longer depends 
on q. Note also that we do not make any assumption on 
the explicit form of 6(t), which, for a given time enters 
Eq. (|2|) only as parameter. 

The lines shown in Fig.lHlare a fit of Eq. Q to the ex- 
perimental data taken for three different values of t with 
A^o, b and qo as free parameters. As one can see our sim- 
ple model reproduces the experimental points very well 
within the error bars which are taken as the standard 
deviation of a five measurement average. In Fig. ^ we 
report the values of qo obtained from the fits shown in 
Fig. ini together with the theoretical prediction obtained 
for a homogeneous cylindrical condensate. As expected, 
the value of qo measured for the different times does not 
exhibit a significant dependence on the BEC-lattice inter- 
action time t (see Fig.EJi). We can directly compare these 
values with the theoretical prediction for the threshold of 
energetic instability for a homogeneous cylindrical con- 
densate I23. Given our uncertainty on the density of 
the sample, which propagates into the theoretical calcu- 
lations, the agreement between theory and experiment is 
good. Furthermore one can extract an estimate for the 
characteristic time of energetic instability by looking at 
the decay of the number of atoms above threshold b{t). 
Assuming for b an exponential decay (i.e. b = bo e~^^^^^) 
the timescale for energetic instability for the experimen- 
tal condition of Fig. |31 is tei ~ 400 ms, as shown in 
Fig. Due to the experimental difficulties in control- 
ling the thermal fraction of the initial atomic cloud, we 
did not study the dependence of tei on the condensed 
fraction. This would be a very interesting measurement 
since it could provide a further insight into the role of 
the thermal fraction. However, as we pointed out in the 
description of Fig. 13 as the condensed fraction increases 
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FIG. 3: Number of atoms remaining in the condensate as a 
function of quasimomentum in a lattice with s=0.2 for differ- 
ent interaction times t. The three datasets refer respectively 
to t = 150 ms (circles), t = 300 ms (triangles) and t = 600 ms 
(squares). The initial cloud has a condensed fraction of about 
65%. The curves are fit to experimental points using Eq. @ 
with A^o, b and as free parameters. 



the number of atoms remaining in the condensate after 
a fixed t increases as well. 



this dissipation using a radiofrequency shield in order to 
make the measurement with no discernible thermal frac- 
tion. This is particularly important if we want to sepa- 
rately address the two regimes of energetic and dynam- 
ical instability. The results of this procedure are shown 
in Fig. O where we plot the BEG lifetime as a function of 
quasimomentum for s = 0.2 with and without RF-shield: 
without the use of an RF-shield (open circles) it is im- 
possible to distinguish the onset of dynamical instability 
8it q = 0.56 qb because this feature is masked by the 
reduction of the number of atoms caused by energetic in- 
stability. On the other hand, by keeping an almost pure 
condensate throughout the experiment (filled diamonds 
in Fig.ISI) the discontinuity in the lifetime entering the dy- 
namically unstable regime becomes clearly visible despite 
a residual reduction of the lifetime due to non complete 
efficiency of the RF-shield. An important point is that, 
even though without RF-shield there is no measurable 
threshold, the lifetime deep in the dynamically unsta- 
ble regime is the same with or without the RF-shield as 
dynamical instability is the dominant loss mechanism in 
both cases. We will show in the next section, devoted to 
a deep experimental investigation of this regime, that in- 
deed dynamical instability takes place much faster than 
energetic instability. We note here that, as already well 
established in literature l27|, the introduction of an RF- 
shield increases the lifetime even in the absence of the 
optical lattice and this explains the difference in the mea- 
sured lifetime at g = 0. 
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FIG. 4: (a) Fitted values for at different interaction times 
(symbols) and theoretical prediction for a homogeneous cylin- 
drical condensate with the same experimental parameters 
(solid line). Dashed lines mark the error on the theoretical 
value which is due to uncertainty in the measured peak den- 
sity of the sample, (b) Number of atoms above given by the 
product A^o^ obtained from the fit to the experimental data 
as a function of interaction time t: continuous line is a fit to 
an exponential decay from which we measure a characteristic 
time of tei = (416 =b 55) ms. 



These results demonstrate that a thermal fraction trig- 
gers the dissipative mechanism connected with energetic 
instability and therefore it is possible to strongly reduce 
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FIG. 5: Lifetime of the condensate with and without RF- 
shield as a function of quasimomentum for the first Bloch 
band (logarithmic plot). The lattice height is s = 0.2 and 
vertical lines are drawn in correspondence of the calculated 
thresholds for a homogeneous cylindrical condensate respec- 
tively for energetic (dotted) and dynamical (dashed) instabil- 
ity 21]. Error bars are smaller than point size. 



6 



0.05n 




0.5 

Quasimomentum [q/qj 

FIG. 6: Measured loss rate of the BEG (reciprocal of lifetime) 
as a function of quasimomentum in the first Bloch band with 
s = 0.2 and s — 1.15 respectively. The vertical lines corre- 
spond to the theoretical values for the threshold of dynamical 
instability 21]. The horizontal error bar is due to the presence 
of the confining harmonic potential as it will be explainded in 
the Appendix. 



IV. "ZERO TEMPERATURE" 
MEASUREMENTS: DYNAMICAL INSTABILITY 

For q > 0.5 qb we observed the onset of dynamical in- 
stability characterized by two distinct signatures which 
we will discuss in some detail: a threshold value of quasi- 
momentum above which the lifetime dramatically de- 
creases (as pointed out above) and, for higher values of 
quasimomentum, the presence of complex structures in 
the density distribution of the expanded atomic cloud 

In Fig. ini we plot the reciprocal of the lifetime of the 
condensate (loss rate) as a function of quasimomentum 
in the first Bloch band for two different values of lat- 
tice height. For both values of s we observe a precise 
value of quasimomentum for which there is a sudden in- 
crease in the loss rate of atoms from the BEC. For our 
experimental parameters, this threshold is almost inde- 
pendent on the condensate density so that, differently 
from the case of energetic instability, inhomogeneity does 
not play a significant role in this context. Indeed there 
is a remarkable agreement between the experiment and 
the calculated thresholds for a homogeneous cylindrical 
condensate shown in Fig. Elas vertical lines. 

From a theoretical point of view the onset of dynamical 
instability is signalled by the appearance of excitations 
of definite wave vector /c, that grow exponentially in time, 
thus modifying the momentum distribution of the system 
characterized by peaks spaced by the periodicity of the 
lattice in the reciprocal space, 2qB |1,[2^- The character 
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FIG. 7: Calculated growth rates |Im(a;)| of the excitation 
modes of an infinite cylindrical condensate in a periodic lattice 
with s = 1.15, as a function of the condensate and excitation 
quasimomenta q and k respectively (first Bloch band). 



of these excitation modes can be obtained by considering 
small deviations from the condensate wavefunction and 
solving the corresponding Bogoliubov equations, as dis- 
cussed in Ti, ^29] for the ID case and in 21] for a 3D 
system. For the system considered in this work when 
q ^ 0.5 qB, the frequency of some modes in the excita- 
tion spectrum develops a nonzero imaginary part Im(cj) 
(see Fig. [7|, which is the distinctive feature of dynam- 
ical instability and corresponds to the fact that, once 
excited, these modes will grow exponentially in time. As 
it is shown in Fig. [3 for a given value of q there exists a 
range of such unstable wavevectors /c, which is more and 
more extended as q approaches the band edge. Among 
all these unstable modes the one with the highest growth 
rate plays a major role in the dynamics of the unstable 
condensate. 

The unambiguous attribution of the observations to 
the phenomenon of dynamical instability has been possi- 
ble by comparing the measured loss rates (inverse of the 
lifetime) to the calculated growth rates of the most dy- 
namically unstable modes. The results of this comparison 
are shown in Fig. [51 for a lattice with s = 1.15 . We point 
out that a full quantitative comparison between theory 
and experiment cannot be performed because the mea- 
surements occur outside the validity of the linear analysis 
on which the theory is based. However, the distinctive 
agreement between the theory and the experiment in- 
deed shows that the observed phenomenon is dynamical 
instability and that the mode which is most unstable in 
the very initial stage of the dynamics (i.e. when linear 
theory is correct) imprints its timescale on the following 
dynamics. 

For those quasimomenta for which dynamical insta- 
bility is more severe and thus lifetimes are shorter, we 
observed the appearance of complex structures in the ex- 
panded density profiles, shown in Fig. [31 for s = 1.15 
and q = 0.55 ^b- These can be the signature either of a 
density modulation or a phase fragmentation that leads 
to the observed fringes through an interference-like effect 
30] . As already stated in we have observed the dis- 
appearance of the fringes as the condensate reverts to its 
initial state if we let it evolve in the magnetic potential 
alone, after switching off the lattice. This process how- 
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FIG. 8: Comparison between the measured loss rates of atoms 
from the condensate (reciprocal of lifetime) and the theoret- 
ically calculated growth rates of the most unstable modes in 
the linear regime. Lattice height is s — 1.15. Note that the 
measurements span the first three Bloch bands. 



ever takes place in a much longer timescale (~ 800 ms) 
than the one of dynamical instability (lifetime can be 
of the order of a few milliseconds). In order to com- 
pare our observations with the theory, we have simulated 
the actual experimental procedure by solving the time- 
dependent 3D GP equation 



zfi^^(x,t) = 



'2m 



V' + F(x, t)^gN\^\ 



^(x,t), 



(10) 

where V is the sum of the harmonic ^ and periodic 
potential From the solution of Eq. (|10|) we extracted 
the axial power spectrum, defined as (the tilde indicates 
the Fourier transform along z) |20| 



P{p,) = 27r rdr\^{r,p,)\' 



(11) 



which is shown in Fig. ^1 The upper frame (Fig. Unk ) 
represents the momentum distribution at the end of the 
initial ramp, and is characterized by sharp peaks local- 
ized at integer multiples of 2qB as discussed above. Then, 
in accordance with the prediction of the linear analysys. 
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FIG. 9: (a) Absorption images of the expanded condensate 
after different interaction times with a lattice with s = 1.15 
for two different values of quasimomentum. Note the sudden 
change of timescale crossing the threshold of dynamical in- 
stability at Q = 0.525 and the appearance of structures in 
the density profiles for the unstable case {q = 0.55 ^s). (b) 
Reabsorption of excitations following 5 ms of interaction with 
the lattice and different times of evolution in the pure har- 
monic potential after switching off of the lattice. In all these 
pictures the lattice moves from top to bottom. 



some modes of complex frequency start growing as indi- 
cated by the arrow in Fig. ITOb (in the numerical simula- 
tion these modes are triggered by the numerical noise). 
Afterwards, the nonlinear dynamics introduces processes 
of mode mixing and the momentum distribution gets 
complicated fFigs. ITOh.d). However the structure is still 
characterized by well localized peaks close to qb^ and the 
position of the peaks still shows invariance under transla- 
tion of2qB' From these results it is possible to calculate 
the expanded density profile and thus make a direct com- 
parison with the experimental observations. This com- 
parison is shown in Fig. ^3 foi" s = 1.15, q = O.bbqB, 
Ir = 10 ms and t = 50 ms. The simulation well re- 
produces the structure of the central peak observed in 
the experiment, but the naive expectation that the ex- 
panded density profile of the condensate would simply 
reflect the structure of the momentum spectrum is not 
correct. We note that indeed the occurrence of dynami- 
cal instability may lead to a rapid population of the non- 
condensate (thermal) fraction, whose behavior in the lin- 
ear regime is governed exactly by the same Bogoliubov 
equations as for the semiclassical fluctuations of the con- 
densate discussed so far 1^3 • This means that the 
Gross-Pitaevskii approach may fail when dynamical in- 
stability yields a strong depletion of the "coherent" con- 
densed fraction. In order to account for this, one should 
include in the theory also the interaction between the 
condensate and noncondensate fractions, which may be- 
come macroscopically populated for later times . The 
formation of a thermal component could strongly affect 
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FIG. 11: Comparison between the simulated ballistic expan- 
sion of the momentum distribution shown in Fig. llOH (top) 
and the cross-section of the measured expanded density dis- 
tribution of the condensate (bottom). In both cases the ex- 
pansion time is 28 ms, s = 1.15, = 0.55 g^s, t = 50 ms and 
the lattice is adiabatically switched off ramping down its in- 
tensity in 2 ms. Vertical lines mark the three main peaks of 
the calculated momentum distribution corresponding respec- 
tively to 0, 1.01 and 2 KqB- In both cases the lattice is moving 
in the direction of positive z. 
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FIG. 10: Calculated axial momentum distribution P{pz) of 
the condensate in the combined harmonic trap plus optical 
lattice for t = 12,33,48,51 ms (from (a) to (d)) (logarithmic 
plot). The simulation is performed with s = 1.15. The arrow 
in (b) marks the appearance of the unstable modes. 

the expanded density distribution masking the momen- 
tum component populated by the interaction with the 
optical lattice. Actually the experimental density pro- 
file reported in Fig. lTTT bottom) shows an uniform tail on 
the right side of the main peak compatible with a thermal 
fraction dragged by the lattice. 

V. CONCLUSION 

We reported on the experimental observation of ener- 
getic and dynamical instability of a BEG in a moving 
ID optical lattice. A clear separation of the two regimes 
is obtained by controlling the temperature of the system 
and adiabatically loading the condensate in a Bloch state 
with precise quasimomentum. 

On one hand we have shown that energetic instabil- 
ity is deeply connected with the presence of a dissipa- 
tive mechanism such as the one provided by thermal 
atoms around the condensate. We have derived a simple 
phenomenological model to take into account the effects 
of inhomogeneity of the atomic density distribution and 



found a good agreement between this model and the ex- 
periment. On the other hand we have observed that in 
the regime of dynamical instability the lifetime of the 
condensate is critically dependent on quasimomentum. 
We have compared this distinctive dependence of the con- 
densate loss rate on quasimomentum with the theoretical 
prediction on the growth rate of unstable modes in the 
initial regime of the dynamics. We have found a very nice 
agreement thus demonstrating that the observed phe- 
nomenon is dynamical instability and that the initial ex- 
citations play a dominant role in the following evolution. 
Furthermore we compared the structure observed in the 
expanded density profile with the results of the solution 
of time-dependent 3D Gross-Pitaevskii equation finding 
a significant agreement between observation and theo- 
retical calculation. We have been able to make a direct 
comparison between the timescales of the mechanisms of 
instability: we have experimentally demonstrated that 
in the regime of low lattice height and for a mostly con- 
densed cloud, dynamical instability is one order of mag- 
nitude faster than energetic instability. Finally we note 
that the critical effect of a thermal cloud on the lifetime 
of the condensate in the presence of energetic instability 
suggests that it is possible to use a moving optical lat- 
tice to detect the presence of a non-condensed fraction 
well beyond the sensitivity of the imaging system usually 
employed in this field. 
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APPENDIX: THE EFFECT OF THE HARMONIC 
POTENTIAL 

As explained in section ^ the adiabatic loading of the 
condensate in a moving optical lattice can be viewed as a 
slow transformation of the spectrum of the system from 
the free-particle parabola to a Bloch band. During this 
procedure the wavepacket changes from that of a free- 
particle to a Bloch state characterized by a velocity which 
is related to the energy spectrum by the well known re- 
lation Vn,q = dqEn,q^ whcrc En^q is the energy of the 
eigenstate with quasimomentum q in band n. As it is 
shown in Fig. E| the difference between q and the free 
particle velocity increases with increasing q. Since the 

0.8 I < \ < \ < \ < y\ < 1 



where z is the direction of the lattice, m is the mass of a 
Rb atom, qo is the initial quasimomentum at which the 
condensate is loaded and F is the external force acting 
on the atoms. In our case this is the harmonic restoring 
force F = —mu'^z. For low values of where the band 
has a parabolic shape (i.e. En^q = fi^g^/2m*), equations 
()A.1|) can be analytically solved leading to 

z(t) = -xf^^ smixf^u^t] 
Q{t) = go + ^ (l-cos(^c.,t)) 

where vq = hqo{l/m* — 1/m). This solution is an 
oscillation both in real and quasimomentum space at 
the trap frequency rescaled by the effective mass m* = 
^{dqVn,q)~^- The amplitude of this oscillation in real 
space (~ 1/im) is too small to be measured by our imag- 
ing system while the amplitude of oscillation in momen- 
tum space cannot be neglected. Increasing when the 
parabolic approximation of the band fails Eqs. (lA.lj) can 
be still integrated giving asymmetric oscillations. Even- 
tually, as shown in Fig. for sufficiently high values 
of q^ the solution of Eqs. (|A.1|) ceases to be oscilla- 
tory: the condensate motion becomes unbounded and 
quasimomentum indefinitely grows. We took into ac- 




Quasimomentum [q/q ] 



FIG. 12: Plot of the Bloch velocity as a function of quasimo- 
mentum in the first band for s = 1.15. The dashed line is the 
free particle velocity (i.e. v = Tiq/m). Notice that the devia- 
tion of the Bloch velocity from the free particle one increases 
as q increases. 



qo=0.3 qg qo=0.6 




t [s] t [s] 

FIG. 13: Time evolution of a Bloch wavepacket calculated 
from semiclassical Eq. (|A.1|) for s = 1.15 and two different 
quasimomenta qo = O.SqB and qo = O.GqB- 



loading procedure conserves this means that the con- 
densate acquires a finite velocity in the laboratory frame 
and starts moving in the harmonic potential. The center 
of mass motion in the frame of the laboratory follows the 
semiclassical laws of motion for a Bloch wavepacket in an 
external force field [3^ 

f = Vn,q - VL = Vn,q - ^qo/m . ^ 

yhq = F \ ' ) 



count the motion of the Bloch wavepacket in quasimo- 
mentum space considering an uncertainty on q for data 
reported in Figs. El and |H1 given by the amplitude of os- 
cillation in q space, as obtained by numerical solution of 
Eq. (|A.1|) . These considerations explain also why we did 
not take data for values of q close to the band edges for 
which dynamics due to the harmonic confinement was 
not bounded. 
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